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Spectral Function Analysis on Spin Dynamics in 
Double-Exchange Systems with Randomness 
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Spin excitation spectrum is studied in the double-exchange model with randomness. Applying 
the spin wave approximation and the spectral function analysis, we examine excitation energy 
and linewidth using analytical as well as numerical methods. For small wave number q ~ 0, 
the excitation energy is cosine-like and the linewidth shows a q-linear behavior. This indicates 
that the spin excitation becomes incoherent or localized near q — 0. Crossover takes place 
to marginally-coherent regime where both the excitation energy and the linewidth are propor- 
tional to q 2 . The incoherence is due to local fluctuations of the kinetic energy of electrons. 
Comparison with experimental results in colossal magnetoresistance manganites suggests that 
spatially-correlated or mesoscopic-scale fluctuations are more important in real compounds than 
local or atomic-scale ones. 

KEYWORDS: colossal magnetoresistance manganites, double-exchange model, randomness, spin wave excita- 
tion, linewidth 



Introduction — One of main issues in colossal magne- 
toresistance (CMR) manganites AMnC>3 is whether the 
double-exchange (DE) mechanism is enough or not, an 
if not, what is necessary as an additional mechanism 
The metallic ferromagnetism is well described by the DE 
mechanism qualitatively f$ which indicates that the DE 
interaction is obviously an essential element in these ma- 
terials. Recently, however, many experimental aspects 
have been indicated the necessity of additional elements 
beyond it. For instance, insulating behavior above the 
Curie temperature Tc as well as large residual resistivity 
is observed in compounds, such as (La,Ca)Mn03, which 
show relatively low Tea The DE model predicts metal- 
lic conductivity both above and below Tc as observed in 
compounds which have high To, such as (La,Sr)Mn03c) 
Thus, another mechanism might be necessary to describe 
the former 'low-T c ' compounds. Since these deviations 
from the canonical DE behavior appear to be systematic 
for the A-site substitution,!! the question is what is the 
controlling parameter. The magnitude of CMR effects 
becomes larger in lower- Tc materials, therefore the issue 
attracts much attentions from the viewpoint of the basic 
mechanism of CMR. 

One of crucial tests for an additional mechanism is spin 
dynamics. Spin excitation spectrum shows different be- 
havior between high-Tc and low-Tc manganites. In high- 
Tc compounds, the spin excitation shows a cosine-like 
dispersion,E' which is well reproduced by the DE mecha- 
nism alone.EP On the contrary, in low-Tc compounds, 
the spectrum deviates from this form and exhibits 
some anomalies .such as broadening, softening and gap- 
openingooEIJET Several mechanisms h a yp hppn pro- 
posed to explain these anomalies.E3'@B@@@ Most 
of them can reproduce qualitative features of the anoma- 



lous spin excitations to some extent, hence, more quanti- 
tative comparison between experimental and theoretical 
results is desired to determine an essential mechanism. 

An intriguing point is that spin excitations in low- 
Tc manganites show lafgg, hjiJtanaic lincwidths even at 
the lowest temperature.B&EiiO'Lij) This indicates that 
magnons in the ground state are not the eigenstates of 
the system. Theoretical models which can be effectively 
mapped to Heisenberg spin systems with uniform ex- 
change couplings cannot describe this situation. 

The authors have claimed that the anomalous spin ex- 
citation is well reproduced by introducing the random- 
ness.EZP Randomness appears to be promising among 
many proposed scenario since it is inherently controlled 
by the A-site substitution and gives a comprehensive un- 
derstanding for systematic changes of the spin excita- 
tion from the high-Tc to the low-Tc materials, at least, 
in a qualitative level. In the presence of randomness, 
magnons are no longer the eigenstates, and the spin ex- 
citation shows a deviation from the cosine-like dispersion 
and a finite linewidth. It is strongly desired to compare 
the excitation spectrum including the linewidth quanti- 
tatively with experimental results. 

In this Letter, we study the randomness effect on the 
spin wave spectrum quantitatively, especially on the ex- 
citation energy and the linewidth, within the spin wave 
approximation. We find an incoherent magnon excita- 
tion in the vicinity of the zone center and a crossover 
to a marginally-coherent regime where the linewidth is 
almost linearly scaled to the excitation energy. 
Formulation — We investigate spin dynamics in the DE 
model with quenched randomness. The Hamiltonian in- 
cludes both diagonal and off-diagonal disorder, which is 
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explicitly given by 

w=-x;*ii(4 



io^3° 



-h.c.) - Jh E CTi ' Si + E £ i c l 



(1) 

where the first term denotes the electron hopping, the 
second one is the Hund's-rule coupling between itinerant 
electrons and localized spins with a magnitude S, and 
the last one gives the on-site potential. The diagonal 
and off-diagonal disorder are incorporated in the on-site 
potential Si and the transfer integral , respectively. We 
call the former 'the on-site randomness' and the latter 
'the bond randomness' hereafter. In the following, we 
first discuss these randomness effects on spin dynamics 
in a general formulation which does not depend on cither 
details of the type of randomness or the dimension of 
the system. Later, we compare the analytic results with 
numerical ones. 

By applying the spin wave approximation in the lowest 
order of 1/5 expansion, the spin wave excitations in the 
limit of Jn/t — > oo are obtained from the static part of 
the magnon self-energy 
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where f n+ is the fermi distribution function for up-spin 
states. Here tp na {i) is the ra-th orthonormal eigenfunc- 
tion, which satisfies £\ Wjjffe, £i})¥wC?) = (E n - 
aJn)ip na -(i) for Hamiltonian (|l]) with a given configu- 
ration of randomness. The spectral function for the spin 
wave, A(q,bS), is calculated by averaging the quantity 
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for random configurations. Here, N is the system size; 
uii and tpi(j) are the eigenvalues and eigenfunctions of 
Hij, respectively, which satisfy 



(4) 



Note that the observable quantities are to be averaged 
finally for configurations of randomness. 

Before going into the analysis on the spectral func- 
tion, we here discuss the magnon self-energy in de- 
tail. Eq. (0) can be written in the form 
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when we define Bji = ^ n f n t Pn(j) i Pn(^) anc ^ use ^ ne re ~ 

lations Em^m^mWVmO') = J2mk n ^k^m(k)^(j) = 

Hij and f n E n ip n (j)ip^(i) = J2 k KjkBki with the or- 
thonormal property of <p n (i)- Here we drop the spin in- 
dices for simplicity. Apparently from eq. (JsJ) , liij satisfies 
the sum rule 
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Moreover, the matrix element consists of the transfer 



energy of electrons as 
1 



H 



1=5 



25 



, a ^ij ( c i C 3 ) 



25 



25 J %■ 



e ' 



(7) 



Hi 



where Jy = (c|cj)/45 2 is the exchange coupling of the 
corresponding Heisenberg model within the lowest order 
of the spin wave expansion^ The bracket denotes the 
expectation value in the ground state for a given config- 
uration of randomness. Hence the following summations 
equal to the kinetic energy of electrons as 

En ri = -En y - = -(T)/25, (8) 

i i=£j 

where T is the first term in the Hamiltonian (Q). 

Now we apply the spectral function analysis on the 
spin excitation spectrum. By using the m-th moment 
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the excitation energy lu sw and linewidth 7 of the spin 
wave excitation are obtained by 



M = ^q\ 



i\q) = n 
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(11) 



respectively. The moments in eq. (^|) can be calculated 
by the magnon self-energy Hij. By using eqs. (|3|) and 
(0), we obtain 
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Namely, Qg is a Fourier transform of n™. This spectral 
function analysis is valid only if the excitation spectrum 
is single-peaked. Previous study shows that this is the 
case .as long as q ~ even in the presence of random- 
ness.EZP This will be demonstrated also in Fig. 1 later. 
Thus, we can discuss the excitation energy and linewidth 
at q ~ by this spectral function analysis. 

First, we analyze the excitation energy (10). Here- 
after, we assume the electron hopping only for nearest- 
neighbor sites in model (|l]) for simplicity. Then, by 
eq. (j^), Hij has nonzero matrix elements only for i = j 
and j + T). Here 77 is a displacement vector to the nearest 
neighbor site. From eqs. (10) and (|l2"|), the excitation 
energy is given by 
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where we define the site- averaged quantity 11(77) = 
jj Yli^-i,i+T] an0 - use the sum rule (|6|). The symmetry 
for the direction of 77 is expected to be recovered after 
the random average, therefore we assume 11(77) = — A 
to be irrespective of 77. Finally, the spin-wave excitation 
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energy is expressed as 

c sw ( <7 )=A^(l-e i ^). (14) 
V 

Thus, for hypercubic lattices, the spectrum has the co- 
sine form for small q even in the presence of disorder. 
The quantity A gives the spin stiffness since 

LU sw {q)~Aq 2 (15) 

in the limit of q = \q\ — > 0. The spin stiffness is propor- 
tional to the kinetic energy of the electronic Hamiltonian 
(|l|) in our analysis as easily shown by eq. (||). 

Next, we consider the linewidth (11). Similarly to 
eq. (13), the second moment is written as 

a q= E n a (»h,ffe)(l- e i9 ^)(l -e 1 ^), (16) 
V1V2 

where n 2 (»7i, V2) = jf Si ^i+rj^i^i+Vi ■ In the absence 
of disorder, we have Tli+r/.i = ^i,i+rj = — A. Hence the 
linewidth 7 becomes zero since = {n { q } ) 2 . In the 
presence of disorder, the linewidth becomes finite. For 
simplicity, we consider a special q such as (q, 0, 0, • • •) 
on a hypercubic lattice. Denoting two different cases of 
q 1 = r) 2 and r) l = -r] 2 as 

^2(Vi,V 2 ) = A 2 ±^A 2 for77 1 = ±T7 2 , (17) 
we obtain 

= 4(A 2 (l-cosg) 2 + <SA 2 sin 2 <j). (18) 
In the limit of small q, the linewidth is estimated as 

7 2 ( 9 ) ~ 46A 2 q 2 + (A 2 - \sk 2 - A 2 )q\ (19) 

Therefore, there is a q-linear contribution in the 
linewidth 7 in the presence of disorder. 

Consequently, we obtain u> sw cx q 2 and 7 cx q in the 
limit of q — » 0. This indicates that the spin wave exci- 
tation becomes incoherent or localized in the vicinity of 
q = since w sw < 7. This incoherent behavior comes 
from local fluctuations of the transfer energy of electrons 
since the coefficient of the q-linear term in 7 is propor- 
tional to 

SA 2 cx n 2 (j7,77) - n 2 (j7, -77) cx ^ {^i+riA - ^h-rf^) 

i 

00 y] { J i+ri,i - Jj-nA ■ ( 20 ) 

i 

As increasing q, the incoherent regime is taken over by 
the marginally-coherent one where uj sw cx 7 cx q 2 . We 
will discuss this crossover in comparison with numerical 
results in the following. 

Results — For further understanding, we calculate the 
excitation energy and the linewidth numerically, and 
compare them with the above analytic expressions. In 
numerical calculations, we obtain the spectral function 
A(q,u>) following the method in ref. [l?], and calculate 
the moments by the definition (^). In the following, 
we discuss two-dimensional cases as an example. The 
qualitative feature of the spectrum is almost indepen- 
dent of the dimension of the system. The system size 



is 40 x 40 sites and the electron density is fixed at 
n = J2i( c l c j)/N — 0.7. The random average is taken 
for 50 different realizations of random configurations. 
We examine effects of the on-site and the bond random- 
ness separately. In both cases, we consider the binary- 
alloy-type distribution of randomness as £j = ±fe and 
tij = t ± St, where the sign takes plus or minus in equal 
probability. We note that other type distributions, such 
as Gaussian-type, do not change conclusions. As an en- 
ergy unit, we use the half bandwidth W = 1 at the 
ground state for Jh = Se = St = 0. The dimension- 
less constant 2S is set to be unity. 

Figure 1 shows a typical excitation spectrum in the 
presence of randomness. The gray-scale contrast shows 
the intensity of the spectrum, namely, the magnitude of 
the spectral function A(q, lj). The global structure keeps 
a portion of the cosine-like one in the case of pure sys- 
tcm,Lr however, some anomalies appear near the zone 
boundaries- This is due to the Friedel oscillation of Jy 
in eq. (7).EZP The gray curve shows the excitation energy 
^sw(g) calculated from eq. (10). As clearly seen, w sw (q) 
well describes the excitation spectrum in the small-q 
regime where the anomalies are not substantial. This 
strongly supports the applicability of the above spectral 
analysis with eqs. (^)-(ll) in this regime. 

The inset in Fig. 1 shows the numerical results for 
uj S w{q) and -f(q) from q = (0,0) to (n,0). They follow 
u>sw cx q 2 and 7 cx q in the small-q region as predicted 
in eqs. ( |l5| ) and (pj|). At q ~ 0, there is the incoherent 
regime where 7 > o; sw as discussed above. 

In Fig. 2, we compare effects of the on-site and the 

bond randomness by plotting the linewidth as a function 

of the excitation energy. For the bond randomness, the 
1/2 

behavior 7 cx o; S w dominates the spectrum as shown in 

the figure. On the contrary, for the on-site randomness, 
1 11 

the Ol> S w part is observed only in the small-q region, and 
the marginally-coherent behavior 7 cx lo sw is dominant in 
the wide region of q. Thus, the q-linear contribution in 
the linewidth is more dominant in the case of the bond 
randomness compared to the case of the on-site one. This 
indicates that the bond randomness tends to make the 
spin excitation more incoherent than the on-site one. 

1 /2 

As discussed above the 7 cx u> S w behavior is due to lo- 
cal fluctuations. Therefore, this different aspect between 
the on-site and the bond randomness indicates that the 
bond randomness induces larger local fluctuations than 
the on-site one. This is qualitatively understood as fol- 
lows. There are two factors which contribute to the local 
fluctuations of the transfer energy (^J^fji — Ji—Ti,i) in 
eq. (20). One is the expectation value (c\cj) and the 
other is the transfer integral iy. In the case of the on- 
site randomness, we have uniform = t, hence only 
(cjcj) contributes to the local fluctuations. The expec- 
tation value (c\cj) is disordered, however, it may have 
some correlations due to the Friedel oscillation. By the 
Friedel oscillation, the charge density shows a correlation 
whose length scale is ~ 2fcp, where k-p is the Fermi wave 
number. On the contrary, in the case of the bond ran- 
domness, Uj fluctuates from site to site independently in 
model (0). This may suppress the correlation in (cjcj). 
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Therefore, we have stronger local fluctuations and more 
incoherent spin excitation in the case of the bond ran- 
domness than in the on-site one. 

Comparison with experiments — In real materials, some 
spatially-correlated or mesoscopic-scale randomness may 
exist due to, for instances, A-site clustering or twin struc- 
ture of lattices. Our results suggest that if these kinds 
of correlated randomness are substantial, the incoherent 
regime diminishes and the marginally-coherent behavior 
7 oc o; sw dominates the spectrum since they may reduce 
local fluctuations. Experimentally, the incoherent behav- 

1/2 

ior 7 oc u S w has not been clearly observed yet. Perring 
et al. have reported an approximately linear scaling of 
7 oc lo sw in the wide region of the spectrum in layered 
perovskite manganites. B They observed a significantly 
large damping, that is, ^/lo sw ~ 0.3 — 0.5. Our results 
in the marginally-coherent regime in Fig. 2 (a) are con- 
sistent with these experimental observations. Therefore, 
we speculate that the randomness plays an important 
role on the spin dynamics in CMR manganites, and that 
correlated or mesoscopic-scale randomness might be sub- 
stantial compared to atomic-scale one. 

Another possible origin for the linewidth broadening 
which survives down to the. lamest temperature is the 
magnon-electron interaction.ESEj In this scenario, the 
linewidth is predicted to be 7 oc q d+3 where d is the di- 
mension of the system. For d = 2 or 3, this q dependence 
is stronger than in our r esult (pi), and is inconsistent 
with the above experimental result. 
Summary — We have studied the spin dynamics in 
the double-exchange model with quenched randomness 
within the spin wave approximation in the lowest order of 
1/5 expansion. We have derived analytic expressions for 
the excitation energy and the linewidth by the spectral 
function analysis. Incoherent magnon excitation is found 
in the vicinity of the zone center where the linewidth 
is proportional to the square-root of the excitation en- 
ergy. As increasing the wave number, this incoherent 
behavior is taken over by the marginally-coherent one 
in which the linewidth is proportional to the excitation 
energy. Atomic-scale randomness enhances the incoher- 
ence through the local fluctuations of the kinetic energy 
of electrons. Comparison with experimental results is 
satisfactory and suggests an importance of correlated or 
mesoscopic randomness in real materials. 

This work is supported by "a Grant-in- Aid from the 
Mombukagakushou" . 
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Fig. 1. Spin excitation spectrum in the case of the on-site ran- 
domness 5e = 0.25. The gray curve is the excitation energy 
calculated by eq. (10). Inset: The excitation energy oj aw and the 
linewidth 7 as a function of (g, 0) for < q < it. The lines are 
the fits by q 1 and q for u sw and 7, respectively. 
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Fig. 2. The linewidth plotted as a function of the excita- 
tion energy in the case of (a) the on-site randomness Se = 
0.05,0.1,0.15,0.2,0.25 and (b) the bond randomness St = 
0.025,0.05,0.075,0.1,0.125 from bottom to top, respectively. 



